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g . A one-dimensional Ising model with nearest neighbour interactions is applied to study com- 

paction processes in granular media. An equivalent particle-hole picture is introduced, with the 
holes being associated to the domain walls of the Ising model. Trying to mimic the experiments, 
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a series of taps separated by large enough waiting times, for which the system freely relaxes, is 

considered. The free relaxation of the system corresponds to a T = dynamics which can be 

analytically solved. There is an extensive number of metastable states, characterized by all the 

holes being isolated. In the limit of weak tapping, an effective dynamics connecting the metastable 

CN ■ states is obtained. The steady state of this dynamics is analyzed, and the probability distribution 
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function is shown to have the canonical form. Then, the stationary state is described by Edwards 
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thermodynamic granular theory. Spatial correlation functions in the steady state are also studied. 
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I. INTRODUCTION 



Granular systems have attracted the attention of physicists in recent years. A review of 
some features of granular matter can be found in Refs. [1], ^|. One of the most outstanding 
problems is the phenomenon of compaction, i. e., the increase of the density of a loosely 
packed granular system when submitted to vibration or, more generally, to some kind of 
external excitation. Compaction has been extensively analyzed in a series of experiments 
by the Chicago group || [|, Starting from a low-density configuration, near the random 
loose packed state, a system of monodisperse glass beads was vertically tapped. Between 
taps, a long enough time was considered, so that the system reached a mechanically stable 
(metastable) configuration before the next tap started. The density was measured in the 
metastable states, and its evolution as a function of the number of taps was studied. The 
parameter controlling the dynamics of the system is the dimensionless vibration intensity 
T = a/g, where a is the peak acceleration in the tap, and g is the gravity. The density was 
observed to increase very slowly towards a steady value following an inverse logarithmic law 
0, |J, the steady density being a monotonic decreasing function of the vibration intensity 
|5|, 01. Several models, with different underlying physical mechanisms, have been proposed to 
understand this behaviour [0, |[ || [K], [EL], 0, but a complete and detailed physical theory 
is still lacking. 

Simple Ising systems are often used as a first approximation to many different phenomena 
in statistical physics. In a recent work [OL [1^] , the one- dimensional Ising model with nearest 
neighbour interactions has been applied to analyze the problem of compaction in dense 
granular media. This system can be regarded as one of the simplest cases of spin models 
on random graphs, which have been very recently used to investigate different aspects of 



granular matter [|ll], y| [16[ . The "tapping" process is simulated in the following way, in order 
to mimic what is done in the experiments [[3], f|, |5| . First, the system freely relaxes until it gets 
stuck in a metastable configuration (n = 0). This is done by considering a modified Glauber 
dynamics at T = 0, in which only those transitions lowering the energy of the system are 
permitted. Thus, all the configurations with no spins antiparallel to both of their nearest 
neighbours are metastable, i. e., they are absorbent states |DJ of this T = dynamics. 
Starting from the metastable state n = 0, the system is "vibrated" by allowing each spin to 
flip with probability p, independently of the state of its neighbours. Afterwards, another free 
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relaxation at T = is considered, and the system reaches a new metastable configuration 
n=l. By repeating this process, a chain n = 0,1,2,... of metastable configurations is 
generated. It is found that the energy of the system is a monotonic decreasing function of 
the number of taps n. This Ising model can be mapped on a particle-hole model, in which a 
particles is associated to site % if spins % and i + 1 are parallel, and a hole is spins % and i + 1 
are antiparallel. In this way, holes correspond to the domain walls between arrays of parallel 
spins. Interestingly, the dynamics at T = can be analytically solved in the particle-hole 
picture [18[]. As a decrease in the energy corresponds to a decrease in the number of domain 
walls, the density of holes is a monotonic decreasing function of the number of taps n in the 
tapping process, i. e., compaction takes place. In the reminder of the paper, we will refer to 
this system as the TIM (tapped Ising model). 

An analogous description of tapping processes was previously introduced for the one- 
dimensional one-spin facilitated Ising model (1SFM) |10|, [19| , originally proposed by 



Fredrickson and Andersen in the context of structural relaxation in glassy systems [20 



In this model there are no interactions, but only an applied external field h, and a spin can 
only flip if at least one of its nearest neighbours is in the excited state. This system is also 
equivalent to a particle-hole model, in which particles are associated to the spins aligned 
with the field, while holes correspond to the spins in the excited state. Then, facilitated 
spin flips are equivalent to adsorption and desorption of particles on the one-dimensional 
lattice. These processes can only occur provided that there is a hole on at least one of the 
nearest neighbour sites of the flipping spin. The tapping process is modelled in the following 
way |l(|. First, the system relaxes at T = 0, where only adsorption processes are allowed, 
until the system reaches a metastable configuration n = 0. Due to the facilitation rule, the 
metastable states are characterized by all the holes being isolated, i. e., surrounded by two 
particles. Then, starting from the metastable state n = 0, vibration is introduced by let- 
ting the system evolve, with both facilitated adsorption and desorption, for a given time to. 
Afterwards, the system relaxes again at T = 0, which leads the system to a new metastable 
state n = 1. Iteration of this process gives a set n = 0, 1, 2, • • • of metastable configurations, 



with the density of particles being an increasing function of the number of taps n JTU], |T9|j . 

One of the most interesting physical questions in the problem of compaction is the de- 
scription of the steady state reached by the system in the limit of an infinite number of taps. 
Note that thermal energy is irrelevant for granular systems. The important energy scale 
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for a grain of mass m and diameter d is mgd, where g is the gravity. In a typical granular 
system, mgd/kBT ~ 10 12 at room temperature. Therefore, while molecular systems explore 
phase space due to thermal fluctuations, in a powder thermal energy is negligible. Unless 
the system is externally perturbed, each metastable configuration would last indefinitely. 
Thus, thermodynamics is not directly applicable to powders. Nevertheless, some years ago, 
Edwards and coworkers ^] made the hypothesis that the steady state of an externally 



perturbed granular system can be described by an extension of the usual statistical mechan- 
ics concepts to granular media. The central point is the ergodic hypothesis for externally 
perturbed powders: in the steady state, all the metastable configurations of a granular as- 
sembly occupying the same volume are equiprobable. Besides, its most stable configuration 
corresponds to the minimum volume. Therefore, the volume of a granular system is the 
analog to the energy for a molecular system. The entropy is defined as the logarithm of the 
number of metastable configurations, which is expected to be an extensive quantity. Then, 
it is possible to define a new parameter, the compactivity X = dV/dS, playing the role of 
the temperature in a molecular system, with the limit X = giving the most compact state. 

In the last years, a lot of effort have been carried out in order to understand if the above 
"equilibrium statistical mechanics" or "thermodynamic" approach describes accurately the 
steady state of an externally perturbed granular system. Most of this effort has been fo- 
cused on the analysis of simple models [OL |TJ, [TP], ^3], |23|, |27] , although there has 
been also some attempts to test Edwards' theory in experiments with real granular systems 
f|. Very recently, the theory has been checked in a numerical experiment with a realistic 
granular matter model, specially conceived to be reproducible in the laboratory f28|. In the 
context of simple models, there is some numerical evidence of the validity of the thermo- 



dynamic description in the limit of gently tapped systems [13, p4L E5L E6L E7[, although for 



stronger tapping the situation is not clear. In fact, numerical results p4[ |27| show that in 
the limit of strong tapping the Edwards measure does not provide an accurate description 
of the stationary state, at least in some spin models. This has suggested an extension of 



the Edwards approach by introducing a "restricted" measure [[24], £7|. On the other hand, 
analytical results are scarce, even for the simplest models. To the best of our knowledge, 
the only system in which Edwards theory has been analytically derived for tapping dynam- 
ics is the 1SFM described above, in the limit of weak tapping [15 1. An effective dynamics 
for the tapping process, connecting metastable configurations, was obtained and the steady 
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probability distribution was shown to have the canonical form. This leads to a relationship 
between Edwards' compactivity and the dimensionless vibration intensity. Let us also men- 
tion that Crisanti et al. [^9j have studied one- dimensional kinetically restricted models to 
address the validity of the Stillinger- Weber construction, an approach that is related to the 
Edwards measure. 

Due to the lack of analytical results, it seems interesting to investigate the possibility of 
deriving a thermodynamical description in the steady state of other simple models. This is 
an important task from a theoretical point of view. Firstly, it is a relevant question if the 
thermodynamic picture is valid or not for models reproducing the experimentally observed 
behaviour in granular systems. Secondly, if the answer is positive, it might be possible 
to derive relationships between the parameters controlling the evolution of the system, for 
instance the tapping intensity V in the compaction experiment, and the compactivity X, 
which characterizes the stationary state. In this paper, we will center on the analysis of 
the effective dynamics and the steady state of the TIM, also in the limit of a gently tapped 
system. We will obtain the effective dynamics between metastable states as an analytical ap- 
proximation to the original tapping dynamics. This will allow us to derive, also analytically, 
the steady state probability distribution. 

The paper is organized as follows. In Sec. |I| the model is introduced, while the analytical 
solution of the modified Glauber dynamics at T = is presented in Sec. fTJ. Section [TV] 
is devoted to the derivation of the effective dynamics in the limit of weak tapping. The 
properties of the steady state reached in the limit of an infinite number of taps are discussed 
in Sec. |V|. It is shown that the steady probability can be written in the canonical form, with 
the role of the energy played by the volume and the temperature being substituted by a new 
parameter, the compactivity, which is related to the tapping intensity. Then, the steady 



state follows the statistical mechanics theory of Edwards and coworkers |2T], |22[. Finally, 



Sec. [VI| contains a summary of the work and some final remarks. 



II. AN ISING MODEL AT T = (FREE RELAXATION) 

The hamiltonian of the one-dimensional Ising model is given by 

N 



n = -j^2(T i (T i+1 , (2.i) 



1=1 
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where J > is the ferromagnetic coupling constant, N is the number of spins on the lattice, 
and o>i = ±1 is the spin variable at site i. We will consider periodic boundary conditions, so 
that formally <Jn+i — a \- The time evolution of the system is governed by single-spin-flip 
dynamics |JIJ. The probability p(cr,t) for finding the system in configuration cr = {a"j} at 
time t obeys a master equation of the form 

N 

d tP {(T, t)=J2 [w(tr\Ri(r)p(Ritr, t) - w(R i( r\<T)p(<T, t)} . (2.2) 

i=l 

Here Ri<r is the configuration obtained from cr by just changing the state of spin i, and 
w{a\a') stands for the transition rate from configuration cr' to cr. Following Lefevre and 
Dean fl~3| , H| , we introduce a T = dynamics such that only those spin flips decreasing the 



energy of the system are allowed. Namely, the transition rates are 

a 

w(RiCr\cr) = — (1 - o^io*) (1 - 0iC7 m ) ■ (2.3) 

It is easily verified that the above expression vanishes unless spin i is antiparallel to both of its 
nearest neighbours. The constant a defines the basic time scale of the system. Interestingly, 
the most general transition rates bringing the Ising model to equilibrium at temperature T 
are 
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w(Ri(T\cr) 



a 
1 



1 + (5crj_i(J l+ i 



1+5 , . , f 2 J 

— - — Oi \pi-\ + <J i+ i) tanh 



k B T 



(2.4) 

where 5 is an arbitrary constant. The usual Glauber dynamics corresponds to the choice 
5 = 0, whereas Eq. Q2.3p is the zero temperature limit of the case 5 = 1. The dynamics 
defined by the transition rates fl2.3|) cannot be solved in the standard way ||30|1 , since the 
hierarchy of equations for the moments C n = (o"j(jj + i) is not closed. 

Let us go to an equivalent description of the Ising model in terms of particles and holes, 
by introducing a new set of variables 

m = • (2.5) 

If spins at sites % and i + 1 are antiparallel, it is rrii = 1 and we will refer to the site % as 
empty or, equivalently, as being occupied by a hole. On the other hand, if spins at site i 



and i + 1 are parallel it is rrii = 0, and site i is occupied by a particle. Therefore, holes are 
associated to the boundaries between domains of parallel spins, i. e., to the so-called domain 
walls of the system. Due to the periodic boundary conditions, the number of holes must be 
even in any configuration. 

In terms of the variables, the hamiltonian Q2.1p reads 

N N 

H(m) = -J 5^(1 - 2mi) = -JN + 2J^ mi . (2.6) 

i=i i=i 

A dimensionless energy per spin e can be defined as 

-7F--! + Jit <•*>■• M 

i=l 

where we use the notation 

(A(m)) t = J2M™)p(m,t), (2.8) 

m 

for an arbitrary function A(m) of the site variables raj. 

In the hole-particle description of the dynamics, the flip of spin % modifies both the values 
of rrii-i and mj. Then, from Eq. ( |2.3| ) we get for the transition rates in the m variables 

w(Ri^iRim\m) = anii^xnii , (2-9) 

where Ri is now the operator which transforms m,; into 1 — m 8 . The master equation for the 
probability p(m, t) is 

d t p(m,t) = [w(rn\Rj-iRirn)p(Ri- 1 R i rn, t) 

i 

—w(Ri-iRiin\m)p(m,t)} . (2.10) 



In the dynamics defined by Eqs. (|2.9|) and (2.10|), the only possible transitions are the 



simultaneous adsorption of two particles on any two neighbouring empty sites. After a long 
enough time, the system becomes trapped in a metastable state characterized by all the 
holes being isolated, i. e., all the empty sites surrounded by two particles. Of course, the 
particular metastable state reached by the system will depend on the initial conditions. 

The present model displays some similarities as compared with the one dimensional one- 
spin facilitated Ising model (1SFM) [|(J at T = 0. In the latter, a particle can be adsorbed 



on an empty site as long as at least one of its nearest neighbouring sites is empty [10, |19 
Although the dynamics of both models are not equivalent, the metastable states are the 
same, being characterized by having all the empty sites isolated. 



III. ANALYTICAL SOLUTION OF THE DYNAMICS AT T = 

Let us define the set of moments 

D r (t) = (m k m k+1 ■ ■ ■ m k+r ) t , (3.1) 

with r > 0. In the following, we will restrict ourselves to homogeneous states, so that D r (t) 
does not depend on the position k of the first site considered. The lowest moment 

D (t) = (m k ) t (3.2) 

is the density of holes. This quantity can be related to the energy per particle by means of 
Eq. (|3), 

e = -l + 2D . (3.3) 
A hierarchy of equations for the moments D r (t) is easily obtained from the master equation 

d t D r (t) = -arD r (t) - 2aD r+1 (t) , (3.4) 
valid for all r > 0. In order to solve the above hierarchy we introduce the generating function 

oo r 

G(x,t) = Y,^D r (t), (3-5) 

from which all the moments D r (t) can be obtained through 

~d r G(x,t)~ 



D r (t) 



dx r 



(3.6) 



J x=0 

From Eq. ( |3.4j ), it follows that the function G(x,t) obeys the first-order partial differential 
equation 

d t G{x,t) +a(x + 2)d x G(x,t) = 0. (3.7) 
By using standard techniques, the general solution of the above equation is found to be 

G(x,t) = G [(x + 2)e~ at - 2] , (3.8) 

where the function G (y) is the initial condition, i. e., 

oo r 

G Q {y) = G{y^)=Y, V - D ^)- ( 3 - 9 ) 



r=0 



In the long time limit it is 

G(x, oo) = G {-2) , (3.10) 

so that 

lim D {t) = G (-2) (3.11) 

and 

lim D r (t) = (3.12) 

for r > 1. The last result shows that all the holes become isolated in the long time limit, and 
the probability of finding r + 1 consecutive holes, which equals D r , vanishes for r > 1. The 
asymptotic density of holes depends on the initial state, as indicated by Eq. In fact, 



it is trivially seen that the hierarchy of equations (|3.4j ) admits as a solution any constant 
value of Dq as long as D r = for r > 1. 

Now, let us specify the initial condition. We will consider that the system is in equilibrium 
at temperature T at t — 0. The equilibrium distribution is given by the canonical distribution 
corresponding to the hamiltonian (|2.1|). 



or, equivalent ly, 

P ' (m)= E m e- w ' (3 ' 14) 

where (3 = (ksT)' 1 . Using the hamiltonian ( |2.6| ), it is easy to show that 

N 

p e T q (m) = l[p^(m l ), (3.15a) 



e -2(3Jmi 



Pt(jth) = 1 + e _ 2fiJ - (3.15b) 



Therefore, 



-20J 



D e Q \ = (m k ; T) eq = - ; _ 9flJ = a , (3.16a) 



D l q T = (m k m k+1 ■ ■ ■ m k+r ; T) eq = (m k ; T)^ 1 = a r+1 , (3.16b) 

where < a < 1. The value a = corresponds to (3 — > oo (T — > + ), and a = 1 to (3 — > — oo 
(T — > _ ). In the limit (3 — ► (T — > oo), it is a = 1/2, which corresponds to the most 
disordered state. 
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Therefore, the initial condition corresponding to an equilibrium state is given by 

D r {0) = a r+1 , (3.17) 

which leads to 



oo 



r 



G (x) = G(x, 0) = V — a r+1 = ae ax . (3.18) 

r=0 

With this choice, the solution given by Eq. (|3.8|) takes the form 

G(x,t) = aexp{a[(a; + 2)e- at -2]} (3.19) 

and, consequently, by using Eq. (|3.6|), 

D r (t) = a r+1 exp [-art + 2a (e~ at - l)] . (3.20) 

As pointed out above, all the moments D r with r > 1 vanish in the long time limit, while 
the asymptotic value of the density of holes reads 

lim D (t) = ae~ 2a , (3.21) 

i^oo 

which depends on the initial density of holes a, being always smaller than it, since only 
adsorption processes are allowed in the T = dynamics. The dimensionless energy per 
particle in the metastable final configuration follows directly from Eqs. ( |3.3j ) and ( |3.21|) , 



-l + 2ae - ^. (3.22) 



This expression agrees with the result obtained by Lefevre and Dean |13], P|. The energy 
Sao is maximum for a = 1/2, i. e., when the system starts from a fully random configuration. 
Of course, this is equivalent to say that D (oo) has a maximum for a = 1/2. The existence 
of this maximum is in contrast with the result for the 1SFM, where the asymptotic density 
of holes is a monotonic function of the initial density JTD], |3T . 



Therefore, at T = the following picture emerges. Starting from any configuration, 
the system evolves until all the holes become isolated, i. e., it gets stuck in a metastable 
state characterized by all the moments D r vanishing for r > 1. Going back to the spin 
description, the metastable states are those such there is no spin antiparallel to both of its 
nearest neighbours. In other words, all the domains of parallel spins have, at least, a length 
of two sites. The density of holes in the metastable state, or the density of domain walls in 
the spin image, depends on the initial configuration, being given by Eq. ( |3.21| ). 
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IV. EFFECTIVE DYNAMICS FOR TAPPING PROCESSES 



Let us consider the model introduced in section || to get a theoretical approach to the 
compaction processes in vibrated granular systems. The model is tapped in the following 



way |13| , T3| , trying to mimic what is done in the experiments with real granular materials 
H ||, |!|. Firstly, the system freely relaxes as described in the previous section, until getting 
trapped in a metastable configuration (n — 0), characterized by the absence of spins being 
antiparallel to both of their nearest neighbours. These configurations are absorbent states 
|T7| for the dynamics at T = 0. Secondly, starting from the metastable configuration n = 0, 
the system is "vibrated" . Each spin can flip with a probability p, independently of the state 
of its neighbours. Afterwards, another free relaxation at T = follows, until the system 
becomes again stuck in a new metastable configuration, n — 1. By repeating this process, a 
chain of metastable configurations n = 0, 1, 2, ... is generated. It is important to note that 
we can restrict ourselves to values of p in the interval < p < 1/2, since the same evolution 
of the energy is obtained for both p and 1 —p. This is because a probability 1 — p is equivalent 
to a simultaneous flip of all the spins (which does not change the energy), followed by a flip 
of each spin with probability p. Making use of the notation introduced in the introduction, 
we will refer to this tapped one-dimensional Ising model as the TIM. 

If p <C 1, the evolution of the system is very slow, since it is very improbable that a given 
spin flips. The dynamics will be dominated by those transitions in which only a few spins 
change their state during a vibration process. Therefore, for p< 1 an expansion in powers 
of p may be useful, since we hope that retaining the lowest orders would provide a good 
approximation. 

In the previous section we have shown that, in the description of particles and holes, 
the metastable states are characterized by having all the holes isolated, i. e., there are no 
domains of parallel spins with length / < 2. We are going to consider the evolution of the 
system in a single tap, defined as the sequence of a vibration process followed by a free 
relaxation to a metastable configuration, to the lowest order in p. In this limit, only one flip 
during the vibration takes place. The analysis to be presented depends on the length I of 
the domain of parallel spins containing the flipping spin. Our goal is to obtain expressions 
for the transition probabilities W e {(m'\m) from the initial metastable configuration m to 
the final one m' in a single tap. 
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Let us assume first that the flipping spin belongs to a domain of initial length I = 2, i. 
e., the transition is like 

• ••UTTU-..— ►.■4UTU--- , (4.1) 

where the flipping spin has been underlined. The probability that this process occurs in 
the cluster above is p(l — p) 5 . Afterwards, in the free relaxation at T = the up-spin has 
to flip necessarily, reaching a new metastable state in which the domain with I = 2 has 
disappeared. Therefore, to first order in p 

Wrf(uuumnu) = 2p, (4.2) 

where the factor of 2 follows because of the other path connecting the same initial and final 
states, and corresponding to the flip of the spin on the right of the domain with I = 2 
during the vibration. These trajectories are shown in FIG. [I], both in the spin and in the 
particle- hole pictures. Introducing an usual notation, in the transition rates W c i{m'\m) we 
have only indicated the sites involved in the given rearrangement. In the particle and hole 
picture, it is 

Wef(00000|01010) = 2p , (4.3) 

the process involves the elimination of both holes. During the vibration, one hole diffuses 
next to the other one, so that in the free relaxation two particles are simultaneously adsorbed 
on them, as it is shown in FIG. [[]. 

When the flipping spin belongs to a domain of length I > 2, we have to analyze two cases 
separately, by distinguishing whether the flipping spin during the vibration is located at the 
domain wall or not. In the former case, for instance 

l 

...unr*~^---uin---, (4.4) 

no transition occurs in the following free relaxation, since the final state in the vibration 
process is metastable. Therefore, the effective transition probability between the metastable 
states is 

WefUun ium)=p. (4.5) 

Similarly, 

w^fCTTili I TTTIi) =p- (4.6) 
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FIG. 1: (a) Destruction of a domain of length I = 2 in the tapping process. In the vibration, one 
of the spins of the domain flips with probability p, and the remaining unstable domain of length 
/ = 1 disappears in the free relaxation at T = 0. (b) The same trajectories, in the particle and 
hole picture. 



These transitions are one-site diffusion processes of a hole in both directions, 

Wef(0010|0100) = Wef(0100|0010) =p, 



(4.7) 



which are present for alH > 2. 

When one of the internal spin flips, we have to analyze separately I = 3 and I > 4. First, 
let us consider 1 = 3. The internal spin is a nearest neighbour of both the spins at the 
domain walls, and the rearrangement occurring in the vibration has the form 



• ••lit nil... — ...utitu... 



(4.8) 



Afterwards, in the free relaxation, any of the three spins of the domain can flip, with the 
same probability, i. e., 1/3. If it is the central spin the one nipping, returning then to its 
original state, nothing has occurred globally, and the group is again in the initial state. On 
the other hand, if any of the external spins of the domain flips first, say the one on the 
left, the reached state is not metastable yet. Then, the free relaxation involves another flip, 
in which the upward spin has necessarily to go down. Therefore, the effective transition 
probability for the complete process is 



WefCillllll | IITTT1I) = -p. 



(4.9) 



The trajectories leading to this rearrangement are shown in FIG. |[ The factor of 2/3 
appears because the final state is the same, independently of which is the first external 

13 



1L 



(b) 




10010 



p 



* 1 1 1 1 



000000 



FIG. 2: (a) Trajectories leading to the destruction of a domain of length I = 3 in a tapping process. 
In the vibration, the central spin flips with probability p. Afterwards, the system freely relaxes 
at T = 0. Paths leading to a final state identical to the initial one are not shown, (b) The same 
trajectories as in (a), in the particle-hole picture. Note that the flip of one spin corresponds to the 
change of two consecutive sites in the associated particle lattice. 

spin nipping in the free relaxation, and each trajectory in FIG. |] contributes p/3. In the 
particle-hole picture, the process consists of the occupation of two holes separated by two 
particles, with a probability 



In a domain of length I > 4, there are / — 2 internal spins. Suppose that the one next to 
the left wall flips in the vibration, 



To the right of the flipping spin, a stable domain of length I — 2 appears, and the spin to its 
left must move downwards in the free relaxation, i. e., the domain wall moves two sites to 
the right in the whole process. The probability of the transition is 




(4.10) 



■ ■■ U T T TT • ■■ 



..•lit ITT-.. 



(4.11) 



VMUUTTIUTTTT)=P, 



(4.12) 



or, in the particle-hole description, 



Wef(00010|01000) =p. 



(4.13) 



The process corresponds to a two-sites diffusion of the hole to the right. Obviously, it is also 
possible that a hole diffuses two sites to the left, corresponding to the flip of the internal 
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TABLE I: Probabilities of the first order transitions in a single tap, connecting metastable states. 



Process Initial state Final state W e f 

One-site diffusion 0100 0010 p 

0010 0100 p 

Two-sites diffusion 01000 00010 p 

00010 01000 p 

Destruction of a hole pair 01010 00000 2p 

010010 000000 \v 



spin next to the right wall in the vibration, 

W e f(01000|00010) = p. (4.14) 

Finally, if the internal spin which flips in the vibration is not next to any of the domain walls, 
it has to return to its original state in the free relaxation, and there is no global transition 
in the tap. 

The obtained transition probabilities up to first order in p are summarized in Table |. 
The particle-hole description is used, since it is more convenient for the analysis of the 
compaction process. Therefore, to the lowest order, with only one flip during the vibration, 
the density of particles cannot decrease and compaction takes place. When higher orders 
are retained, processes leading to a decrease of the number of particles show up, as it will 
be discussed below. 

The description of the effective dynamics between metastable states is not complete, 
even qualitatively, if it is restricted to the lowest order in p. In particular, the existence of 
a steady state characterized by the flipping probability in the vibration p [13, LI] is lost. 



The Markov process describing the dynamics between metastable states is not irreducible, 
and the configuration with all the sites being occupied by particles is an absorbent state 
[T7| of the dynamics. Therefore, in order to have a more complete description, we are led 
to consider higher orders in p, i. e., processes involving more than one transition during the 
vibration. This will be done in a physical way, similar to that of Ref. O. We are not going 



to consider those second order processes (two flips in the vibration) whose effect can be 
obtained by means of a combination of two processes of order p, but only those processes for 
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(a) 

P\ | | 

T T T T T T T T T T 

(b) 

00000^— 10 10 

FIG. 3: (a) Creation of a domain of length I = 2 in the tapping process. In the vibration, two 
consecutive spins belonging to a domain of length I > 6 flip with probability p 2 . As a consequence, 
a new stable domain of length I = 2 shows up. (b) The preceding trajectory, in the particle-hole 
picture. 

which the effective transition probability W e i vanishes to the lowest order. In particular, this 
is the case for all those trajectories decreasing the density of particles, as already mentioned. 
The inclusion of these processes modifies essentially the physics of the tapping process, so 
that they must be taken into account in our effective dynamics. 

Then, we will consider that there are two flips take place in the vibration. As in the free 
relaxation only transitions decreasing the energy are allowed, we have to analyze two cases: 
(i) flip of two nearest neighbours spins, and (ii) flip of two spins separated by one site. If the 
two flipping spins are separated by more than one site, the local free relaxations associated 
to each of them are independent, and the result is a product of two first-order transitions. 

Suppose the transition during the vibration in the cluster depicted in FIG. [3[ whose 
probability is p 2 (l — p) 4 . In the free relaxation no transition can happen, since a metastable 
domain of length I = 2 has been created. Therefore, to the lowest order it is 

^ef(umiiuuu)=p 2 , (4.i5) 

or, in the particle-hole picture, 

iy ef (01010|00000) = p 2 . (4.16) 

In order to derive this transition probability it has been assumed that the domain which 
the nipping spins belong to initially is at least of length I = 6. If it has a smaller length, it 
is easily shown that the resulting rearrangement can be obtained as a combination of two 
first-order processes. 

The other kind of second order processes we have to include in our approximation corre- 
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(a) 



t t t t t t t t t 

,2 



1/3, 

I I * T 



(b) r 1/3 

1111 0^- 10 10 



FIG. 4: (a) In the vibration, two spins of the domain of length I > 7 flip with probability p 2 , 
creating three consecutive unstable domains of length 1 = 1. If, in the subsequent free relaxation, 
it is the central spin the one which flips, a new stable domain of length / = 3 appears, (b) The 
same trajectory, shown for particles and holes. 



sponds to a transition during the vibration of the form 



■ UtltU 



(4.17) 



whose probability, to leading order, is p 2 . In the free relaxation process, first either one 
of the underlined spins which have just nipped or the one between them flips, the three 
changes having the same probability. In the former case, the group of spins returns to the 
initial configuration, all of them upwards. On the other hand, if the central spin flips, a new 
domain of length I = 3 appears, with probability 

1 



Wef(umuiiuuu) = ^ 2 . 



(4.18) 



The trajectory leading to this transition is shown in FIG. f|. In terms of particles and holes 

1 o 



W/ e f(010010|000000) 



(4.19) 



It is easy to convince oneself that there are no more second order transitions increasing the 
number of holes, and that any other second order process can be decomposed in a product 
of first-order transitions. The new transitions appearing to second order are shown in Table 
U Together with the first-order processes in Table |, they define the approximate effective 
model for tapping processes we are going to analyze in the following. It is important to note 
that the introduction of the second order processes is fundamental from a physical point of 
view. The processes described in FIGS. § and|], with probabilities given by Eqs. (|4.16p and 
( |4.19| ), are the inverse of the processes in FIGS. [I] and||, Eqs. (|4.3| ) and ( |4.10|) , respectively. 
Therefore, the Markov process defined by the transition probabilities in Tables | and |I| is 



irreducible [17], i. e., all the states are connected by a chain of transitions with non-zero 
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TABLE II: Probabilities of the second order processes between metastable states leading to an 
increase of the number of holes. 



Process Initial state Final state W e f 

Creation of a hole pair 00000 01010 p 2 

000000 010010 \p 2 

probability. This property will be fundamental in the analysis of the steady state, presented 
in the next section. 



V. STEADY STATE SOLUTION 

As pointed out above, the Markov process defined by the effective transition probabilities 
connecting metastable states is irreducible |HJ and, consequently, there is a unique steady 
state for each given value of p. This steady state will be reached by the system from any 
initial configuration. Besides, it will be shown that the system described by the effective 
master equation verifies detailed balance. By using this property we will be able to obtain 
the steady distribution analytically. With regards to the original model, the expression holds 
in the limit of gently tapped systems, p <C 1, for which the effective transition probabilities 
of Tables | and [Tl] have been obtained. 

In order to calculate the steady distribution function, we will bet a priori on a stationary 
solution P s (m) of the master equation for the tapping process verifying the detailed balance 
condition, 

W ei {m\m')P s {m') = W ci (m'\m)P s (m) . (5.1) 

Given the uniqueness of the steady state, if a solution is found in this way, its own existence 
will be the proof of the detailed balance property in the system. Detailed balance implies 
that all the configurations having the same number of holes k will be equiprobable, 
since they are connected through diffusion processes, which are isotropic. Their probability 
will be denoted by P s (m^). Moreover, for the processes changing the density of the system 
it is 

p g ( m /(fc+2)) ^ W c t(m'V e+ V\mW) _ p 
P s (mW) ~ ty ef (m( fc )|m'( fc + 2 )) ~ 2 ' 
This expression applies for both pairs of transitions with non-zero probability, given by 
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Eqs. flPfl) and (0), and Eqs. Q419|) and ( CT ), respectively (see also Tables | and 0). 
Consequently, 



Z V2, 

where Z is a normalization constant, and we have taken into account that the number of 
holes is always even. Defining a new variable X by 

-i/x = . fp 



\/|, (5-4) 

the steady probability distribution can be written in the "canonical" form 

P s (m«) = V fc /*, (5.5) 

so that X is identified as the compactivity of Edwards' statistical mechanics theory of 
powders JH], 22| . The number of holes plays the role of the volume or, more precisely, the 



excess volume from the densest state. The normalization constant Z is the analog to the 
partition function. From it, all the steady properties of the system can be obtained in the 
standard way. 

The calculation of Z is quite an easy task, 

N/2 

Z = ^2 Z k> (5.6a) 

k=0, k even 



Zk = tt{ N) e- k ' x , (5.6b) 

with being the number of metastable states with k holes for a lattice with N sites. 
The maximum number of holes is N/ 2 (we are assuming that N is even) , and the number 
of holes k must be even in the TIM because of the periodic boundary conditions. A simple 
combinatorial argument leads to 

iZfc " k\(N — 2k)\ ■ 1 ' 

In the large N limit, the sum in Eq. (|5.6a|) can be evaluated by the saddle point method, 

since Zy. has a sharp maximum as a function of k, with the result 

1 i + (i + 4 e -i/xu/2 
In C = ^ In Z = In i±ii±| L_ . (5 . 8 ) 
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0.25 




FIG. 5: Comparison of the numerical evaluation of the steady density of holes, as given by 



Eq. (5.11), (circles) and the analytical explicit expression obtained from the effective dynamics 



approach, Eq. (5.9). The agreement is quite good for p < 0.15. 



The number of holes is the property analogous to the energy of a molecular system, and the 
steady hole density reads 

p,_ dln( (l + 4e-^)V 2 -l 
d(l/X) 2(l + 4e- 1 /*)i/2 ' 1 ' 

It must be stressed that Dq = k/N, being k the value of the number of holes k for which Zk 
reaches its maximum. The steady probability distribution is a very sharply peaked function 
around k, which assures the equivalence of the microcanonical and canonical ensembles 
for the calculation of the mean values of the physical properties in the steady state. The 
stationary density of holes p s = 1— Dq is a monotonic decreasing function of the compactivity 
X. As the compactivity, given by Eq. (|5.4|) , increases with the vibration intensity p, p s is 
also a monotonic decreasing function of the vibration intensity, a behaviour analogous to 
that of real granular systems [H, |J. In the limit p< 1, Eq. (|5.9| ) reduces to 

D s ~ e ~ 1 / X = (5.10) 



In Refs. [L3|, [L4| it was found that the steady density of holes for arbitrary p is given by 
the solution of the equation 

Di = [2 M + Di(l - 4 M) ] exp 

with q = 1 — p. Of course, as Eq. (|5.11| ) gives the exact value of the density of holes, it is 
symmetric against the change p — > 1 — p. In the limit p <C 1, the above expression yields 



(5.11) 
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FIG. 6: Comparison of the numerical evaluation of the steady density of holes, as given by the 
exact solution Eq. ( |5.11| ) (circles), the analytical expression obtained from the effective dynamics 



approach, Eq. ( |5.9| ) (solid line), and their leading behaviour (5.10) in the limit p< 1 (dashed line), 



in the interval 10 6 < p < 10 . It is observed that Eq. ( 5.10| ) only gives a good approximation in 



the very weak tapping regime, p < 10 3 , for which the steady density of particles p s = 1— Dq > 0.98. 



Dq ~ a/p/2, which agrees with the small p limit of the expression obtained from the effective 
dynamics, Eq. (|5.9|) . In FIG. |5], we compare the density of holes as a function of p, obtained 
from Eqs. ( |5.9|) and ( |5.11 ). It is observed that the agreement is quite good for p < 0.15, 



i. e., for a steady density of particles p s = 1 — -Dq > 0.85. On the other hand, the leading 
behaviour for j) < 1, Eq. (|5.10|) , only holds for very small values of p, p ^ 10~ 3 , as it is 
clearly shown in FIG. || Therefore, the accuracy of the results obtained from the effective 
dynamics picture extends further than what might be expected from a second-order theory 
in the nipping probability p. In fact, this is not so surprising, because a similar behaviour 
was found in the one- dimensional facilitated Ising model submitted to tapping processes 



T9| . The wide range of applicability presented by Eq. ( |5.9| ) can be understood on physical 
grounds, by realizing that the effective dynamics approach is not a standard second-order 
expansion in the flipping probability p. Among all the second-order processes, only those 
introducing new physically relevant transitions, i. e., transitions that cannot be written as 
a combination of first-order processes, are collected. This physically motivated expansion 
allows an analytical approach to the rather difficult problem of tapping. 
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The entropy is defined in the usual way, 



S = -J2 p s(m)\nP s (m) = ^D s + \nZ, (5.12) 



X' 

m 

which is an extensive quantity. The specific entropy per site is 

S 1 



N x" 

l + 4e -iA)i/2_i 1 + (1 + 4e-V^)V2 



2X(l + 4e-V^)i/2 



+ In 



(5.13) 



It is possible to define a function analogous to Helmholtz's specific free energy, 

$ = -X\n( = D s -Xa, (5.14) 
so that Eq. ( |5.9|) can be expressed 

■ d($/X) 

D "=wm- (5 ' 15) 

The description of the steady state of the TIM and that of the tapped 1SFM presented in 
Ref. [19] are closely related. In both models, the metastable states are the same, namely those 
characterized by having all the holes isolated, i. e., surrounded by two particles. Moreover, 
in the weak tapping limit, the steady state probability distribution has the canonical form 
in both models. This implies that their "thermodynamical" properties are the same, when 
expressed in terms of the compactivity. Probably, this equivalence does not hold for stronger 
tapping, for which the simple description for the steady state developed here seems to need 



some refinements, as it follows from the numerical experiments reported in Refs. ^4], |27. 
Finally, let us note that the number of holes must be even in the TIM, while the 1SFM is 
free from this restriction. Of course, this difference becomes irrelevant in the thermodynamic 
limit, in which the density of holes is a continuous variable. 



Spatial Correlations 

In the steady state, the only spatial correlations present in the model are due to the 
impossibility of having two nearest neighbour holes. This property can be used to simplify 
the calculation of the correlation functions in the steady state. 
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We introduce two "entities" , a and f3, being a the set formed by a hole with the particle 
that is at its right nearest neighbour site, and (3 a single particle. Thus, any metastable 
configuration is obtained as an, unrestricted, arbitrary permutation of the entities a and (3. 
Here we take a very large system, so that we do not need to consider periodic boundary 
conditions. We will denote by N a and Np the number of entities a and (3 in a given 
configuration, respectively. Therefore, 

2N a + N p = N. (5.16) 

It is obvious that N a equals the number of holes k in the configuration, so that the density 
of entities a and (3 are related to the density of holes D , 

Xa =V^t = D , (5.17a) 

Xf3 = ^ = 1 - 2x a = 1 - 2D , (5.17b) 

where the averages are done over the considered ensemble of systems. 

Since the positions of entities a and (3 are independent in the steady state, it is very 
easy to compute stationary correlations. For instance, the probability of finding two holes 
separated by r particles is given by 

F r = - m k+1 ) • • • (1 - m k+r )m k+r+1 ) s = < r -i Q . (5.18) 

Here, x s a g r -i a is the steady density of clusters composed by two entities a separated by r — 1 
entities (3. The number of entities f3 is r — 1 because the particle in site k + 1 together with 
the hole in site k constitute the first entity a. We will consider that r > 1, since the holes 
are isolated in the metastable states. It is 

x a pr-i a = ^ Na(3 N ^ , (5.19) 

where N a pr-i a is the number of clusters of the kind indicated above. By definition, 

(N a(3 r-i a ) = (N a )P(f3 r - 1 a\a), (5.20) 

where P((3 r ~ 1 a\a) is the conditional probability of finding a cluster composed of r — 1 
consecutive entities (3 and one entity a to the right of one entity a. As the entities a and f3 
are independent in the steady state, the stationary value of this conditional probability is 

P s ([3 r - l a\a) = PsiP^a) = [P^' 1 P * , (5.21) 
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where and PI are the probabilities of finding an a and a (3 entity in the steady state, 
respectively Obviously, 

p.- W - go f5 2 2a) 

p ~ (N a + N )>- l-D-r [ } 

where we have made use of Eq. (|5.17[ ). Therefore, putting together Eqs. ( |5.18 )-( |5".22| ), we 

get 

(D°) 2 fl-2D°\ r . . 

n ' °» <r > 1) . (5.23) 



1 - 2D S Q \ 1 - D s 

Since the moment F r s equals the probability of finding two holes separated by r particles, it 
is clear that 

F o ' = 0, (5.24) 

reflecting that two holes must always be separated by at least one particle. The moments 
F r obey the following "sum rule" , 



J2 F r = D o, (5-25) 



r=0 

expressing that the sum of the probabilities of finding two holes separated by an arbitrary 
number of particles equals the probability of finding one hole, i. e., 

Do = (m k ) = (m k m k+1 ) +(m k (l - m k+1 )} 
= F + Fx + (m k (l - m k+1 )(l - m k+2 )) 

oo 

= ... = ^F r . (5.26) 

r=0 

The calculation of other spatial correlations in the steady state is straightforward, by 
following a line of reasoning similar to the one used to find F*. For instance, F£ also 
provides the probability of finding a cluster composed by two entities a and r — 1 entities 
(3, no matter the way they are ordered, because of the independence of the entities a and (3 
in the steady state. 



VI. FINAL REMARKS 



In this paper we have analyzed a one- dimensional Ising model with nearest neighbour 
interactions formulated in a way appropriated for the study of compaction in granular media. 
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An equivalent particle-hole description has been introduced, in which the holes are associated 
to the domain walls of the original Ising system. The free relaxation of the system is modelled 
by a T = dynamics |13[], which only allows those spin flips decreasing the energy of the 
system. Any configuration with all the holes (domain walls) being isolated is metastable, i. e, 
it does not evolve with this T = dynamics. The tapping process is described as composed 
of two steps: (a) vibration, i. e., starting from a metastable configuration, each spin of 
the system is flipped with probability p, and (b) the system freely relaxes with the T = 
dynamics until it reaches a, in general, different metastable configuration. The parameter 
characterizing the tapping process is the "vibration intensity" p. 

In the particle-hole description, the dynamics at T = is analytically solvable, by writing 
a closed hierarchy of equations for the probability distribution functions D r of finding r + 1 
consecutive holes in the system In the long time limit, the system gets stuck in a state 
where D r = for all r > 1, i. e., all the holes are isolated, as indicated above. 

Tapping is a rather complex process, since each tap is composed of two neatly different 
processes: vibration and free relaxation. In order to get a physical insight into the mech- 
anisms responsible for the behaviour of the system under tapping, the derivation of the 
effective transition probabilities for the Markov process connecting the metastable states 
reached by the system in two consecutive taps is needed. In general, this is a formidable 
task, but in the limit of a gently tapped system these transition rates between metastable 
states can be computed up to the second order in p. 

In the first order, the only possible transitions are diffusion and destruction of a hole pair. 
Then, compaction takes place, since there are no processes decreasing the density of particles 
in the system to the lowest order. To describe the steady state, second-order processes 
must be taken into account, so as to have transitions that increase the number of holes. 
Interestingly, these transitions are just the inverse of those decreasing the number of holes 
to the lowest order. Therefore, the Markov process is irreducible, i. e., all the metastable 
configurations are connected through a chain of transitions with non-zero probability. As 
a consequence, there is an unique, well-defined, steady probability distribution for each 
value of p. Besides, the effective transition rates verify detailed balance. This property has 
been used to derive the steady distribution analytically, finding that it has the canonical 
form. Thus, a relationship between Edwards' compactivity and the vibration intensity p is 
obtained, in the limit of weak tapping. 
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The system analyzed in this paper, as formulated for modelling tapped granular media, is 
closely related to the one-dimensional facilitated Ising model JT0|, [19], In the respective 
particle-hole pictures, the metastable states are the same, those having all the holes isolated. 
Although, in the limit of weak tapping, the corresponding effective dynamics connecting 
metastable configurations are not equivalent, the steady state is described in both cases by 
the canonical distribution. The role of the energy is played by the number of holes and that 
of the temperature by Edwards' compactivity, which is related to the vibration intensity by 
an Arrhenius-like expression. In this way, we find a sort of "minimal" model for Edwards' 
description of the steady state of externally perturbed granular media: a one-dimensional 
system of particles and holes, with the metastable states characterized by having all the 
holes isolated, and a canonical probability distribution function. Nevertheless, in order 
to have an actually complete description of the steady state, a relationship between the 
parameters characterizing the statics and the dynamics of the system, i. e., between the 
compactivity and the vibration intensity, is needed. Then, it is also necessary to derive the 
effective dynamics from the underlying original models (TIM, 1SFM, etc.), when trying to 
understand the steady state behaviour. 

The effective dynamics approach between metastable states has been shown to be a 
powerful tool, in order to study the steady state of models for granular systems submitted to 
tapping processes. It allows to identify the physical mechanisms responsible for the increase 
of the density, and also for the existence of a steady state characterized by a density being a 
monotonic decreasing function of the vibration intensity. In simple models, the calculations 
can be thoroughly done in the limit of gentle tapping, deriving analytically the steady state 
distribution. The results so obtained are consistent with recent extensive numerical tests of 



Edwards hypothesis in simple systems fl3| , [24] , Pq , f2q , p7[| , although the systematic deviations 



found for stronger tapping [24, 27 1 cannot be accounted for within the second-order theory 



developed in this paper and in Ref. [19]. This would need an extension of the effective 
dynamics approach to the whole range of vibration intensities, which is certainly not an 
easy task. Nevertheless, it is hard to believe that the simple structure of the transition 
probabilities of the effective dynamics should remain unaltered for stronger tapping, i. e., 
deviations from the simple canonical distribution found here are to be expected. 
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